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Abstract 

I study the Schrodinger operator with the strong magnetic field, 
considering links between geometry of magnetic field, classical and 
quantum dynamics associated with operator and spectral asymptotics. 
In particular, I will discuss the role of short periodic trajectories. 

Preface 

1 will consider Magnetic Schrodinger operator 

(1) W = ^ (J2 Pjg^'{^)Pk - V) , Pj = hDj - f,Aj 

where g-i^^ Aj^ V are smooth real- valued functions of x G M*^ and {g-'^) is a 
positive-definite matrix, < /? <C 1 is a Planck parameter and // 1 is a 
coupling parameter. I assume that /-/ is a self-adjoint operator. 

2-dimensional magnetic Schrodinger is very different from 3-dimensional, 
all others could be close to one of these cases but are more complicated. 

I am interested in the geometry of magnetic field, classical and quantum 
dynamics associated with operator (1) and spectral asymptotics 

(2) J e{x,x,0)^{x)dx 
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as /? ^ +0, fi +OC where e(x,y, r) is the Schwartz kernel of the spectral 
projector of H and ip{x) is cut-off function. Everything is assumed to be C^. 

1 Geometry of Magnetic Field 

Magnetic Intensity Magnetic field is described by a form 

(3) ^ = ^kdxk) ^j^^'^'j ^ 

k J,k 

with 

(4) Fj, = djV,-dkVj. 

So a does not change after gauge transformation A A + Vcj) and this 
transformation does not affect other objects I am interesting in as well. 

I am discussing local things and Aharonov-Bohm effect which demon- 
strates that knowledge of a, gj^^ V is not sufficient to characterize spectral 
properties of H is beyond my analysis. 

Canonical forms If a is of maximum rank 2r = 2[c//2j one can reduce it 
locally to the Darhoux canonical form 

(5) ^ = XI A dx2j. 

l<i<r 

So, (5) is a canonical form of a near generic point for generic A. However 
situation becomes much more complicated near general point for generic A. 
Complete results are not known. Assuming c/ = 2r and a is generic J. Mar- 
tinet [Ma] had shown that Z/^ = {x, rankF(x) < d — k} are submanifolds 
and calculated their codimensions. In particular, codimZi = 1. Moreover, 
Z2 = as c/ = 2, 4 (not true for d > 6). 

As c/ = 2 generic 2-form a has a local canonical form 

(6) a = xidxi A dx2, Z = {xi = 0}. 

However, as c/ = 4 not all points of Z = Zi are equal: A = {x G Z, Ker F(x) C 
7';^Z} is submanifold of dimension 1. As x G Z \ A dim (Ker F(x) n Tj^Z) = 1 
and in its vicinity one can reduce a to 



(7) 



a — xidxi A dx2 + dx^ A dx^ 
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while in the vicinity of x G A canonical form is 

(8) a = dxi A dx2 — X/^dxi A dx^ + x^dxi A dx/\. + X3C/X2 A dx^ + X4C/X2 A C/X4+ 

2(>^i - ^(x| + X4'))c/x3 AC/X4 

(R. Roussarie [Ron], modified by X2 1-^ X2 — IX3X4). 

Magnetic lines Magnetic lines are described by 

dx 

(9) — GKerF(x)n7.Z 

where one can skip HTxY^ without changing the substance of the definition. 

As rankF = d (and thus d is even) there are no magnetic lines. As 
rankF = d — 1 (and thus d is odd) through each point passes exactly 1 
magnetic line. 

As c/ = 2 and a is defined by (6) the only magnetic line is {xi = 0}. As 
c/ = 4 and a is defined by (7) magnetic lines are straight lines {xi = 0, X3 = 
const, X4 const}. As c/ = 4 and a is defined by (8) A {xi = X3 = X4 = 0} 
and magnetic lines are helices {xi = 0,X3 = rcos9,x/\. — rsin^,X2 = const — 
r^9/2\ (with r = const), winging around A. 

True geometry From the point of view of operator H simultaneous anal- 
ysis of form a and metrics {g^^) should be crucial, but I am not aware about 
any results. It appears, however, that only eigenvalues ±ifj and eigenspaces 
of matrix (F/) = {"^jg'^Fjk) are really important, and in the case of the 
generic magnetic field they are not very difl&cult to examine. 
As c/ = 2 

(10) fi = F,2/Vg, g = det(^^'^)-^ 
while for c/ = 3 

(11) ^1 = ^( E S^'s'-Fj,F,^Y^' = ( Y: Sj^PF'f" 

j,k,l,m j,k,l,m 

where P = ^'^Zk 1 ^ki is a vector intensity of magnetic field, s^^' is an 
absolutely skew-symmetric tensor with s^^^ = V ^/g' 
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2 Classical Dynamics 

2.1 Constant case 

This case has been weU- known long ago. 

2D case Assume first that g^^ , Fj^ and V are constant. Then with no loss 
of the generality one can assume that g-'^ = 6jk^ skew-symmetric matrix (Fj^) 
is reduced to the canonical form: 

{fj j = l,...,r, k=j + r 
-fu J = r+l,...,2r, k=j-r , 
otherwise 

fj>0 and \/ = 0; moreover, one can select Aj(x) as linear functions. 
Then as c/ = 2, /i > classical particle described by Hamiltonian 

(13) H(x, = ^ (E^'^wfe- - ^AW) (^^ - ^^^w) - ^w) 

moves along cyclotrons which in this case are circles of radius pi — {jifi)'^ ^/2E 
with the angular velocity uji — jifi on energy level {AV(x, ^) — E}. 

3D case As c/ = 3, /i > there are a cyclotron movement along circles 
of radii pi = {fifi)~^ y/2E^ with the angular velocity Ui — jifi and a free 
movement along magnetic lines (which are straight lines along Ker F) with 
a speed ^/2Ef and energy E is split into two constant arbitrary parts E — 
El + Ef. 

Multidimensional case Multidimensional case with c/ = 2r = rankF is 

a combination of 2D cases: there are r cyclotron movements with angular 
velocities uoi^ = jif^ and radii = {pfk)~^V2Ek where energy E is split into 
r constant arbitrary parts E — Ei + E2 + • • • + E^- The exact nature of the 
trajectories depends on the comeasur ability of /i, ... , f^. 

As d > 2r = rank E in addition to the cyclotronic movements described 
above appears a free movement along any constant direction v G Ker E with 
a speed ^/2Ef where energy E is split into r + 1 constant arbitrary parts 
E = El + E2 + • • • + E, + E^ 

This difference between cases d — 2r — rankE and d > 2r — rankE will 
be traced through the whole paper. 
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2.2 Full rank case 

Assume now only that d — rankF (see f.e [Ivr2, Ivr3]. In addition, assume 
temporarily that Fj^ and g^^ are constant hut potential V{x) is linear. Then 
cyclotronic movement (s) is combined with the magnetic drift described by 
equation 

(14) ^-^^{2^^)-'Y.<^'^,V 

k 

where = {Fjk)~^- 

As c/ = 2 it wiU be movement along cycloid and multidimensional move- 
ment will be combination of those. 

Not assuming anymore that V is linear we get a bit more complicated 
picture: 

• Equation (14) holds modulo 0{ii~^)] modulo error 0{iLi~^t)] 

• As c/ = 2 cycloid is replaced by a more complicated curve drifting along 
V = const and thus cyclotron radius p = {fifi)~^y^2E + V would be pre- 
served; 

• In higher dimensions all cyclotron radii are preserved as well. 

Without assumption that g-^'^ annd Fj^ are constant picture becomes even 
more complicated: 

• As c/ = 2 cycloid is replaced by a more complicated curve drifting along 
f~^{V + 2E) = const (thus preserving angular momentum uJipl according 
to equation 

(15) f^={2fi)-\\/f-\V + 2E))^ 

where ^ means clockwise rotation by 7r/2 assuming that at point in question 

• In higher dimensions (at least as non-resonance conditions fj ^ Vy ^ k 
and fj ^ fk + fi yy , k, I are fulfilled) one can split potential V = Vi + • • • + V/c 
so that similar equations hold in each eigenspace of {FD and both separate 
energies and angular momenta are (almost) preserved. 
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2.3 3D case 

As c/ > 2r = rankF the free movement is the main source of the spatial 
displacement and the most interesting case is 2r = d — 1 and especially 
c/ = 3, r= 1. 

In this case the magnetic angular momentum M is (almost) preserved; 
thus kinetic energy of magnetic rotation is ^f~^M^] therefore in the coor- 
dinate system such that g"^^ = Sy the free movement is described by ID 
Hamiltonian 

(16) H,{x,,Ci;x',M)^^ei-\v,ff 

with effective potential \/efF(xi,x^) = V — f~^M^^ x — (xi,x^). 

Thus the particle does not necessarily run the whole magnetic line and 
the helix winging around it does not necessarily have constant the step or 
radius. Effect of the magnetic drift is rather minor. 



2.4 2D case: variable rank 

Situation becomes really complicated for variable rankF. I am going to 
consider only c/ = 2, 4 and a generic magnetic form a. Let me start from the 
model Hamiltonian as c/ = 2: 

(17) H' = \(il + {^2-^^xl[/vf-\)■ 

the drift equation is 

and for |xi| ^ ^ = gives a proper description of the picture. 

For the model Hamiltonian (17) with /x = 1 (otherwise one can scale 
xi 1-^ /x^/^xi, ^2 ^ l^k) we can consider also 1-dimensional movement along 
xi with potential 

(19) V(xi; /c) = 1 - (/c - x{/v)\ k = 6; 



Then for odd u 



• V is one-well potential; 
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• As /c = ±1 one of its extremes is where ^(0) = 0; 

• Well is more to the right/left from as ±/c > 0; as /c = well becomes 
symmetric. 



(a) /c> 1; (b) < /c < 1; (c) k = 1; 

Figure 1: Effective potential for odd u. 

On the other hand, for even u potential is always symmetric and 

• We have two- well potential with the central bump above surface if /c > 1 

• and below it as < /c < 1: 

• As /c = ±1 one of its extremes is where ^(0) = 0; 

• Well is more to the right/left from as ±/c > 0; as /c = well becomes 
symmetric. 




(a) /c> 1; (b) /c = 1; (c) < /c < 1; (d) -1< /c < 0; 

Figure 2: Effective potential for even u. 
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Let us consider trajectories on the energy level 0. From the analysis of 
the potential one can conclude that for k ±1 the movement along xi is 
periodic with the period 

(20) 



Jx-{k) ^y2V[Xl■, 



k) 

however one needs to analyze the increment of X2 during this period 

(21) ,(,)=2r*'"(^«i. 

Jx-{k) \/2V{xi\ k) 

One can prove that /(/c) ^ as /c ^ /c* with < /c* < 1 for even v and /c* = 
for odd V. In particular, /c* ^ 0.65 for v — 2. Further, /(/c) x (/c — /c*) as 

On figures 3-5 are shown trajectories on (xi, X2)-plane plotted by Maple in 
the outer zone (these trajectories have mirror-symmetric or central-symmetric 
for even or odd v in zone xi < 0), in the inner zone for even v and in the 
inner zone for odd v respectively. For the spectral asymptotics periodic tra- 
jectories are very important, especially the short ones. Periodic trajectories 
shown above are very unstable and taking V = 1 — axi instead of xi breaks 
them down (figure 6). One can find details in section 1, [Ivr5]. 



2.5 4D case: variable rank 

The most natural model operator corresponding to the canonical form (7) 
is + H" with as above and 2H" = ^| + (^4 - x^f. Then H" is a 
movement integral. Therefore the dynamics is split into dynamics in (x^ — 
{^i^^i^ Ci' C2) described above with potential 1/1/ = V — 2E and the standard 
cyclotron movement with energy E in {x'\^") = (X3, X4, ^3, ^4). 

Situation actually is way more complicated: considering + {1 + axi)E 
we arrive to the 1-D potential V — 1(1 + axi) and playing with E and a one 
can kill the drift even for /c 1 leading to many periodic trajectories. 

Consider canonical form (8) which in polar coordinates in (X3, X4) becomes 

(22) a = c/((xi - \p^)dx2 + (xi - -/)p''de)) . 

The most natural classical Hamiltonian corresponding to this form is 

(23) 2H = il + (6 - ^(xi - \p'')f + + r-\{} - p{x, - ^p^)p^)' - 1 
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3.3-1 
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(a) /c 1; as xi > trajectory 
moves up and rotates clockwise 
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(b) k decreases, still k > 1. Trajec- 
tory becomes less tight; actual size of 
cyclotrons increases; 




(c) k further decreases, still /c > 1; (d) /c = 1. Trajectory contains just 

one cyclotron. 

Figure 3: Movements in the outer zone 



with g, dual to p,9. 

Note that ^2 and d are movement integrals and therefore xi — is 
preserved modulo 0{ii~^). Based on this one can prove that 

• There is a cyclotronic movement with the angular velocity x in the 
normal direction to paraboUoid {— xi + = 
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(a) /c < 1 slightly 



(b) k further decays but still larger 
than /c*. Drift slows down 
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1 / 




(c) k = k* . No drift; trajectory be- (d) k < k"" . Drift now is down! 
comes periodic 




(e) k decays further. Drift down ac- (f) and further; as /c = — 1 we have 
celerates. just straight line down 



Figure 4: Even movements in the inner zone 
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(a) /c < 1 slightly. Drift is up and (b) k decays further, (c) /c = 0. No drift, 
the fastest Drift up slows down Trajectory is periodic 

further 

Figure 5: Odd u] movements in the inner zone 





(a) Even u (b) Odd u 

Figure 6: Breaking periodic trajectories by a linear potential 



• combined in the zone {|xi| < cp^} with the movement similar to one 
described in 2D case in (p, 0)-coordinates (with {xi = 0} now equivalent 
to {p = p}) on the surface of this ellipsoid 

• and also combined some movement along X2; 

• I did not consider zone {|xi| > cp^} since it was not needed for the 
spectral asymptotics. 



One can find details in section 1, [Ivr8]. 
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3 Quantum Dynamics 

Microlocal canonical form {Birghoff normal form) play a crucial role in the 
analysis of the quantum dynamics and spectral asymptotics. 

3.1 Canonical forms. I 

In the case d — 2 and a full-rank magnetic field canonical form of Magnetic 
Schroding operator is (| of) 

(24) uiixi, ^i-HDi){h^Di + i/xl) - l4/(xi, ii-HDi)+ 

m+/c+/>2 

with 6jy = ^- o v|/, 1/1/ = V o v|/ with some map U^. The first line is main part 
of the canonical form. 

In the case c/ = 3 and a maximal-rank magnetic field microlocal canonical 
form of Magnetic Schroding operator is (|) of 

(25) uji{xi^ X2. 11-^ hD2){h'Dl + 11^x1) + 11^01- l/l/(xi, X2, /x-^/7D2)+ 

m+n+/c+/>2 

^2— 2/77— 2/c— /— n^/+A7 

Again, the first line is main part of the canonical form. 

In the case c/ > 4 and a constant rank magnetic field microlocal canonical 
form of Magnetic Schroding operator is of the similar type provided we can 
avoid some obstacles: 

If fj have constant multiplicities (say, fj are simple for simplicity) then 
the main part is 

(26) Yl x"' f^-'hD"){h'Dl^^j + f^'xl^^j) + h'D"- 

l<J<r 

W{x\x',fi-^hD'); 

where x' = (xi, ... , x^), x'' = (x^+i, ... , x^+^), 2r = rank F, q = c/ - 2r. 

Next terms appear if one can avoid higher order resonances: Pjfji^) — 
with Pj G Z; 3 < Ylj \Pj\ is calleed the order of the resonance. 
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After operator is reduced to the canonical form one can decompose func- 
tions as 

(27) u(x) = ""(^'' ^")Tp,(x.+,+i) • • • T,Xxd) 

where T are eigenfunctions of Harmonic oscillator h^D^ + /x^x^ (i.e. scaled 
Hermite functions). 

Then as 2r = c/ one gets a family of r-dimensional /x~-^/7-PDOs and for 
2r < d one gets a family of q-dimensional Schrodinger operators with poten- 
tials which are r-dimensional /x~-^/7-PDOs. 

The similar approach also works for 2 and 4-dimensional Schrodinger op- 
erators with the degenerate magnetic field of the types I considered before but 
only in the far outer zone {7(x) =^ |xi| ^ and to this form operator 

is reduced in balls 6(x, |7(x)). 

3.2 Canonical forms. II 

As c/ = 2, 4 and magnetic field degenerates there is also a more global canon- 
ical form. 

As c/ = 2 in zone {|xi| <C 1} this form is (after multiplication by some 
non- vanishing function) 

(28) h^Dl + (/7D2 - iJix^/iyf - W{x) + perturbation 

with W = (if oiriginally h - (/)dist(x, Z = {xi = 0}). 

For d — A one can separate a cyclotron part corresponding to the non- 
vanishing eigenvalue ^; after this one gets a 3-dimensional second-order DO 
(+ perturbation) with the principal part which is the quadratic form of rank 
2 and a free term V — {2a + l)iihf2 where a G Z+ is a corresponding magnetic 
quantum number. 

This operator could be reduced to the form similar to (28) (at least away 
from A = {x3 = X4 = 0}); here W = {V -{2a + l)iihf2)(^-^'^ . 

Close to A but as |xi| < Cp^ one can get a similar form but with 9 instead 
of X2 and jip instead of ji. 

3.3 Periodic orbits 

One can prove that semiclassical quantum dynamics follows the classical one 
long enough to recover sharp remainder estimates but the notion of periodic 
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orbit should be adjusted to reflect logarithmic uncertainty principle 
(29) |osc(x)| . |osc(D)| > Ch\ log h\ 

where h is eflPective Plank constant (it could be h or fi~^h or one of them 
scaled depending on the particular situation). 

I need a logarithm because I am interesting in the size of the box outside 
of which function is negligible rather than in the mean quadratic deviation. 
Function exp{—\x\^ /2h) scaled shows shows that boxing requires a logarithmic 
factor. 

So instead of individual trajectories I consider their beams satisfying log- 
arithmic uncertainty principle. One can see that the classical trajectory is 
not periodic but cannot say this about the semiclassical beam until much 
larger time, after few "periods" . 




(a) (b) 
Figure 7: Classical (a) and semiclassical (b) periodicity. 



4 Spectral asymptotics 
4.1 Tauberian method 

Problem I am looking at asymptotics as /? ^ +0, /x +oc of 
(30) r((?e)(0) = J (e(., 0)(?;)^^/y = Tr((?E(0)) 
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where e(x,y,r) is the Schwartz kernel of the spectral projector E(r) of op- 
erator H and is a pseudo-differential operator, means a dual operator. 

As Q = / we get TrE(O) which is the number of negative eigenvalues of 
H (and +oc if there is an essential spectrum of H below 0). 

I hope to construct this expession (30) with Q = / from itself for elements 
of the partition of unity with self- adjoint Q > 0. 

Tauberian method Tauberian method Fourier says that the main part of 
r(eQ^) is given by expression 

(31) h-' f (F,^,-i,XT(t)r(u(?;)) dr 

J —oo 

while the remainder does not exceed Cy + C where 

(32) M = Mr = sup | (F,^,-^,XT{t)V{uQl)) \ 

\r\<e ^ ^ 

and s is large, C does not depend on e, T , h, ji and s while C depends on 
e > 0, 7, s. 

Here and below ^(x, y, t) is the Schwartz kernel of the propagator e'^ 
X ^ Co^([-l,l]) equal 1 at [-|, |], x ^ Co^([-l, 1]) equal at [-|, |], 
Xrii) — x{i/T)^ T > and F^^^-i^ is /7-Fourier transform. 

Actually this remainder estimate persists if one replaces T by any larger 
number T' only in expression (31). 

So, I want to increase T without (significantly) increasing Mj in (32). 



Evil of periodic trajectories Microlocal analysis says that if there are 
no periodic trajectories with periods in [y, T\ on energy levels in [— 2e, 26] 
then 

(33) sup I (F,^,-r,XT{t)V{uQl)) I < Ch\ 

|r|<e ^ ^ 

Therefore if there are no periodic trajectories with periods in [7, T'] on 
energy levels in [—26, 26], then one can retain T in (31), Mj in (32) but the 
remainder estimate improves to C-^, + C . 

So, periodic trajectories are one of the main obstacles in getting a good 
remainder estimate. For example, if all trajectories are periodic with the 
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period 7 = 7n then it can happen that Mr x TT^^^Mt^ as 7 > 7n and 
increasing T does not bring any improvement. 

For example, let /x < 1, \/ x 1. Then there are no periodic trajectories 
with periods in [Tq, Tq = Ch\ log/?! and 7"i = e because dist(x(t), x(0)) x 
7 as 7 < 7i and this distance is observable as 7 > 7o. 

Then M < Ch~^To — Ch^~^\ log/?! and the remainder estimate is 
0{h^~^\ log /?!). Actually one can get rid off this log factor in /W = Ch^~^ and 
the remainder estimate becomes 0{h^~^). This remainder estimate cannot 
be improved without geometric assumptions of the global nature. 

Also, taking 7 really small in (31) allows us to calculate u and (31) by 
a crude successive approximation method with an unperturbed operator H 
having coefficients frozen as x = y. 

As /i > 1 the same arguments are true but 7i = e/x~-^ and the remainder 
estimate is 0{iih^~^). This remainder estimate cannot be improved as c/ = 2, 
g"^^, fi and V are constant. 

In our arguments', because all trajectories are periodic (pure cyclotronic 
movement). 

From direct calculations: as domain is all eigenvalues are Landau levels 
{a + |)/i/? — \ V oi infinite multiplicity {a G Z+) and 

(34) e(x, X, r) = 5] e(2r + \/ - (2a? + l)fihf) fih-'f^ 

n>0 

with jumps X fih~^ at Landau levels. 

However in many cases one can improve remainder estimate 0{fih^~^). 
The idea is to show that actually periodicity is broken. 

From the point of view of applications one should take Q with support 
(with respect to x) in ball 6(0, |) (then rescaling arguments could be applied) 
and impose condition on operator only in the circle of light 6(0, 1) with the 
self-adjointness being the only condition outside of it. 

So anything out of 6(0, 1) is a dark territory and we must take 7 < 7* 
which is the time for which trajectory which started from supp Q remains in 
6(0, 1). But we can chose the time direction and we can chose it for every 
beam individually. 

Now, as c/ = 3 the typical trajectory is non-periodic because of the free 
movement and 

• one must take 7 < 7* x 1; 
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• but for most of the trajectories one can take 7"i x 7"* retaining Tq = 



and the remainder estimate is (/?"'" ^) (under very mild assumptions). 
Let d — 2 and fi do not vanish. Then 

• since the drift speed is 0{ii~^) one can take T ^ ji] under certain non- 
degeneracy conditions breaking periodicity of the cyclotronic movement one 
can take Ti — 7* retaining Tq — Ch\ log/?! 

and the remainder estimate is 0{ii~^h^~^). 

When tamed, our worst enemy (periodic trajectories) could become our 
best friend! 

4.2 Results: Constant-rank case 
Results: "Constant" case 

Theorem 1. Let g^^, Fj^ and V be constant and domain be M*^. Then 



where Q.^ is a volume of unit ball in M^. 

In particular, the spectrum is pure point iff r = 0. 

One can prove this theorem easily by direct calculations. Magnetic Weyl 
Expression cE^^(x, E) becomes our candidate for the main part of asymp- 
totics in the general case. 

Results: c/ = 2 As c/ = 2 formula (35) provides a good approximation and 
the non-degeneracy condition below breaks periodicity and provides a good 
remainder estimate [Ivr7]: 

Theorem 2. Let d — 2 and g^^ , Fj^, V be smooth in 6(0, 1); fi non-vanishing 
there and ip G Q^(6(0, 1)). Let assume that all critical values of V/fi are 
non- degenerate. Then 



Ch\\ogh 



(35) ^,^^(x, E) Q,_2.(27r)-^+>^/7- 



X 




(36) 
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as jji < ch . 

Remark 3. (i) If in the general case magnetic field spoils remainder estimate 
0{iih^~^)^ but in the covered case magnetic field improves it; 

(ii) Estimate (36) holds in multidimensional full-rank case as well but non- 
degeneracy condition is pretty complicated and is not generic [Ivr3]; 

(iii) As ei < fih < c nondegeneracy condition changes; as c/ = 2 it reads: 
{2a + l)/x/7 is not a degenerate critical value oi V/f for any a G Z+; 

(iv) As /i + ••• + /",> e > 0, e(x, x, 0) is negligible and = for /x > ch'^. 

Results: Non full-rank case Let us assume that rankF is constant but 
less than d. In this case remainder estimate cannot be better than 0{h^~^) 
but it also cannot be much worse [Ivr4] : 

Theorem 4. Let g^^ , Fjk, V be smooth in S(0, 1), rank(F) = 2r, < 2r < d 
so fi,. . . ,fr do not vanish there and G Q^(6(0, 1)). Then 

(i) As either 2r — d — 1 and some very mild non- degeneracy condition is 
fulfilled or2r = d-2 or ji < h^'^ with (5 > 

(37) I j (e(x, X, 0) - ^,^^(x, 0))^(x) dx\ < Ch'-'; 

(ii) As 2r = d — 1 the left-hand expression does not exceed Cfih^~^~' + Ch'~' 
with arbitrarily small S > 0. 

4.3 Results: Degenerating 2D case 

Consider case c/ = 2, /i x |xi|^~-^ with u > 2 assuming that 

(38) V>eo>0. 

We consider e- vicinity of {xi = 0} with small enough constant e > 0. 

Then in the outer zone {7 = Cfi-^/" < |xi| < e} there is a drift with the 
speed /i~-^7~^, this drift breaks periodicity and therefore contribution of the 
strip {|xi| X 7} with 7 G (7, e) to the remainder estimate does not exceed 
Ch'~' X 7 X ii~'^~^ where the second factor is the width of the strip and the 
third one is the inverse "control time". Then the total contribution of the 
outer zone to the remainder estimate does not exceed the same expression 
as 7 = 7 which is Cfi~'^^h'~'] this is our best shot. 
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In the inner zone {|xi| < 7} or equivalently {|^2| ^ Q} the similar 
arguments work as long as p x — k*V^^^\ > e. 

Otherwise there seems to be no drift to save the day. However it is not 
that bad. Really, period in xi is x 7 and if 

(39) \i2-k*V'/'\xp 

the speed of the drift is x p, the shift with respect to X2 is x p7 and in order 
to be observable it must satisfy logarithmic uncertainty principle P7 x p > 
Ch\ log h\ because characteristic scale in ^2 is p now. So, periodicity is broken 
provided 

(40) P>Pi = C{T'h\\ogh\)K 
which leaves us with much smaller periodic zone 

(41) Z,,, = {\^2-k*V'/^\<p,}. 

And in this periodic zone picking up x 7 we can derive remainder estimate 
0{h~^pi) which does not exceed our dream estimate Cp~^^^h~^ provided 
pi < 7 or 

(42) /i< C(/7|log/7|)-i 

Actually for the general operator rather than the model one we need to 
assume that p > C7 but this does not spoil our dream estimate. 

So we need to consider periodic zone defined by (41) assuming that (42) 
does not hold. 

Inside of Periodic Zone Even in the periodic zone Zper periodicity of 
trajectories can be broken as 1/1/ = V(()~'^/^^^^^xi=o is "variable enough" which 
leads us [Ivr5, Ivr6] to 

Theorem 5. Let c/ = 2, /i = (j){x)\xi\''~-^ with v >2 and condition {38) he 
fulfilled. Then as ^ e CS^{B{0, 1) n {|xi| < e} 

(43) I J (e(x, X, 0) - 5^^^(x, 0))7/;(x) dx\ < C{p-^'' + h^^^^^'^)h'-' 

where here and below h — p^^^h, q = in the general case and q — 1 under 
assumption does not have degenerate critical points '\ 
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To improve this remainder estimate one should take in account the short 
periodic trajectories. ActuaUy, periodicity of the trajectories close to them 
is broken but only after time Tq = Cp~^h\ \ogh\ (see our discussion in sub- 
section 3.3). Skipping details [Ivr5, Ivr6] 

Theorem 6. Let d — 2, fi — ^(x)|xi|'^~-^ with v >2 and condition {38) he 
fulfilled. Then as ^ e CS°{B{0, 1) n {|xi| < e} 



(44) I |(e(x.x,0)-£,^^(x,0))^(x)c/x-|C^(x2,0)^(0,X2)c/x2|< 

provided either some very mild nondegeneracy condition is fulfilled or ji < 
ch^~^ where S^'!^{x2,0) is defined below. 

Remark 7. (i) In the case we are considering right now (and no other case 
considered in this article) condition /"i + • • • + > cq fails and therefore 
e(x, X, 0) is not negligible as /x > ch~^] 

(ii) On the other hand, e(x, x, 0) is negligible as /x > c/?"^; 

(in) As ch~^ < M < c/?""", f'^^(x, 0) is supported in the strip {|xi| < 

7i == Co(/x/7)~-^/^^~-^^} where 7i > 7 = c/x~-^/^; therefore the main part of 
the spectral asymptotics (after integration) is of magnitude {iJ.h)~^/^^~^^ . 

In theorem 6 the correction term is defined by 

(45) £^'^{X2, t) = (27r/7)-^ j no(r; X2, 6, h) c/6 - j S^ir; xi, X2, h) dx^ 
where no is an eigenvalue counting function for an auxiliary ID-operator 

(46) ao(x2, 6- ^) = ^ (r?Dl + (6 - xl/vf - l/l/(x2)) 

and 8^"^ is Magnetic Weyl approximation for the related 2-dimensional op- 
erator. 

Using Bohr-Sommerfeld approximation one can calculate eigenvalues of 
30(^7^27^) with 0{h^) precision and 8^{x2,t) with 0{h~^h^) precision. In 
particular, modulo 0{h~^h) — 0(7"-^) 



(47) £:r(x2,o) = >ch-^n^w\-i^G[^^^) 



1 I 1 
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with some constants x and So and function G defined by 

(48) G{t) = l{t+\v'-[t+ + \\)dv 

with the converging integral in the right-hand expression. One can prove 
easily that 

(49) G ^0, G{t + 1) = G{t), [ G{t)dt = 0,GeC'2. 

Jo 

This is one of examples of the short periodic trajectories really contributing 
to the asymptotics. 

4.4 Results: Degenerating 4D case 

4D case is way more complicated than 2D one. But there are some good 
news: since fi + f2>em the generic case, we need to consider only /x < ch~^. 

The main difficulty in 4D case comes from the outer zone 
1^1 1 ^ because there could be short periodic trajectories. In other words: 
Landau level {2ai + l)fihfi + (2^2 + l)l^hf2 could be flat 0. Still in contrast to 
the very general case when this can happen for up to (/x/?)"-^ pairs a G 
in the assumptions of theorem below it can happen only for no more than C 
pairs and I was able to prove [Ivr8] 

Theorem 8. Let F is of Martinet-Roussarie type and condition {38) be ful- 
filled. Then as is supported in 6(0, 1) n {|xi| < e} 

(50) I J (e(x, X, 0) - S^'^ix, 0))^{x) dx\ < Cii'^'^h-^ + C^^h'^. 

One can improve this result under extra condition breaking flat Landau 
levels [Ivr8]: 

Theorem 9. In frames of above theorem assume that (\///^2)xi=o does not 
have degenerate critical points. Then as is supported in 6(0, l)n{|xi| < e} 

(51) I j (e(x, X, 0) - S^'^ix, 0) - C"!"!^. 0))Hx) dx\ < C/i-^/^/^-s 

Here £^ — 0{is'/^h~^/^) is associated with periodic zone {|xi| < c/i~-^/^}^ 
and is the sum of similar expressions in 2D case for \/^ — \/ — (2/3 + V)\ihf2 
with f3 G Z'^ ; locally all of them but one could be dropped. 
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4.5 Results: Nondegenerating 4D case revisited 

As I mentioned, even if magnetic field does not degenerate, non-degeneracy 
condition of [Ivr3] is not generic (it defines open but not everywhere dense 
set). Recently I was able to prove [Ivr9] 

Theorem 10. Let us consider 4D Schrddinger operator with non- degenerate 
generic magnetic field. Then 

(i) For generic potential V asymptotics holds woth the principal part as in 
(51) and the remainder estimate 0{ii~^h~^); 

(ii) For general potential V asymptotics holds woth the principal part as in 
{50) and the remainder estimate 0{ii~^h~^ + ij}h~^); correction term could 
he skipped with no penalty unless /7~^/^+^ < M < h~^/^~^ in which case it can 
he skipped with the penalty 0{h~^^^~^ where S > is arhitrarily small 

(added February 2, 2008) 

All results here are proven in the series of articles at 
http : // www . math . toronto . edu/ ivr ii/Research/Prepr ints . php 
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